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ABSTRACT 

We analyse the evolution of turbulence and gravitational instability of a galactic disc 
in a quasi-steady state governed by cosmological inflow. We focus on the possibility 
that the coupling between the in-streaming gas and the disc is maximal, e.g., via dense 
clumps, and ask whether the streams could be the driver of turbulence in an unstable 
disc with a Toomre parameter Q ~ 1. Our fiducial model assumes an efficiency of 
~ 0.5 per dynamical time for the decay of turbulence energy, and ~ 0.02 for each of 
the processes that deplete the disc gas, i.e., star formation, outflow, and inflow within 
the disc into a central bulge. In this case, the in-streaming drives a ratio of turbulent 
to rotation velocity a/V ~ 0.2 — 0.3, which at z ~ 2 induces an instability with 
Q ~ 1, both as observed. However, in conflict with observations, this model predicts 
that a/V remains constant with time, independent of the cosmological accretion rate, 
because mass and turbulence have the same external source. Such strongly coupled 
cosmological inflow tends to stabilize the disc at low z, with Q ~ a few, which may be 
consistent with observations. The instability could instead be maintained for longer, 
with a properly declining a/V, if it is self-regulated to oscillations about Q re 1 by 
a duty cycle for disc depletion. However, the 'off' phases of this duty cycle become 
long at low z, which may be hard to reconcile with observations. Alternatively, the 
coupling between the in-streaming gas and the disc may weaken in time, reflecting an 
evolving nature of the accretion. If, instead, that coupling is weak at all times, the 
likely energy source for self-regulated stirring up of the turbulence is the inflow within 
the disc down the potential gradient (studied in a companion paper) . 



Key words: galaxies: evolution 
galaxies: kinematics and dynamics ■ 



- galaxies: formation - galaxies: high-redshift 
galaxies: star formation - methods: analytical 



1 INTRODUCTION 

The basic kinematical properties of galaxies are the aver- 
age rotation velocity V ro t, and the velocity dispersions a 
of their components, namely the random motions of stars 
and the gas turbulence in the interstellar medium (ISM). 
Disc galaxies are supported against gravity by rotation. At 
low redshift, massive discs have gas velocity dispersions of 
~ lOkms" 1 (|Dib et al.l2006T ). with a /V vn t re 0.05-0.1. Their 
stellar velocity dispersion is typically at the level of tens of 
kms -1 and it varies with stellar surface density (jBottemal 
1 19931 : Ivan der Kruitl l20ld : IWestfall et all l201ll ). Observa- 
tions show larger ga s velocity dispersion s at higher redshifts 
(|Epinat et al.l |201C| : iDavies et atll201ll h such that typical 
massive disc galaxies at z re 2 have a re 30 — 80kms _1 and 
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cr/Vr o t re 0.15-0.3 (lErb et alj|2004 iForster Schreiber et all 
l2006l : ICresciet al.ll2009D . 

Galactic discs are often assumed to maintain marginal 
gravitational i nstability, with a Toomre Q parameter Q ~ 1 
(Toomrdll964h . In this case, 



i ;, 



(2v)-^5, 



(1) 



where S is the mass fraction in cold disc within the disc 
radius, 



M dis 
Mtoi 



(2) 



and v is a factor of order unity that depends on the shape 
of the rotation curve, 



1 + 



logKo 

logr 



(3) 
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JDekel et all l2009bh . Eq. JTJ explains why at high red- 
shift the high fraction of cold disc, in terms of gas and 
young stars, requires a high value of cr/Viot in order to 
maintain Q ~ 1. The instability of discs with high 8 and 
c/Viot is more 'violent', in the sense that the structures are 
large r and the dynamical processes are faster (IGenzel et al 



20081 ; iDekel et al.ll2009bl; IGenzel et alJlJOlll ; ICeverino et al 
20121 ; iBournaud et all 1201 lal lbf). The perturbations associ 



ated with the instability, in the form of extended transient 
features and bound clumps, reflect the larger characteristic 
Toomre mass, 



A/clump ^ 1.2 

Mdisc v 



a 

Vrot 



(4) 



Torques between the perturbations drive angular momen- 
tum out and generates mass inflow, partly as clump mi- 
gration and partly as inflow of inter-clump mass. The 
times cale for inflow can be estimated in several different 
ways (|Gammiell200ll ; lDekel et alj|2009bj ) to be 



^dyn 

where 



Vrt 



t,ly 



<r 



R 



(5) 



(6) 



£1 is the angular velocity and R is the effective radius of the 
disc. Thus, a higher a/V I0 t is associated with a faster inflow 
in the disc. 

Since I SM turbulence decays on one or a few dynamical 

1 19981 ; IStone et al.lll998l ; iGammiel 

there must exist a con- 



timescales (IMac Low et al 
l200ll ; lElmegreen fc Scald 
tinuous energy source that maintains the turbulence over 
cosmic time. However, the nature of this energy source is 
highly debatable. The mechanisms that drive the ISM tur- 
bulence could be divided into three kinds. First, stellar 
feedback, such as supernova feedback and radiative feed- 
back from stars, which deposit energy and momentum into 
the ISM. Second, the energy source could be the gravita- 
tional energy released by the instability-driven inflow down 
the potential gradient within the disc, which is a natu- 
ral mechanism for self-regulating the instability at Q ~ 



JWadaet al.ll2002J;lAgertz et al.ll200d; IBournaud et al-lhoid 



iKrumholz fc BurkertJl20ld ; ICacciato et alj201ll ). Third, the 
driver of turbulence could be the kinetic energy transferred 
from the cosmological inflow of clumpy gas. At high redshift, 
the gas streams in as supersonic streams that follow the fila- 
ments of the cosmic web. The streams, consisting of merging 
galaxies and a smoother component, penetrate through the 
halo to the vicinity of the central disc w here they deposit a 
fract i on of their energy and momentum IIBirnboim fc Deke! 



back could be the main driver of turbul e nce in loca l mas- 
sive galaxies (IMac Low fc Klessen 20041; Kim et all l200ll ; 



iDib et all 120061 ; Ide Avillez fc Breitschwerdtl l2007f ) as well 
as in high-redshift galaxies ( Hopkins et al.l 120111 ). On the 



other hand, other estimates and simulations argue that stel- 
lar feedback is unli kely to drive a velo c ity dispersion larger 
than m lOkms" 1 dJoung et all [20091; IDekel et al.l [2009bt 
IBournaud et al.ll2010l ; Idstriker fc Shettvll201ll '). 

Unlike stellar feedback, where the energy emerges 
from nuclear processes inside stars, both the second and 
third mechanisms refer to 'gravitational heating', where 
gravitational potential energy is released as a result 
of infall into a potential well (jDekel fc Birnboiml [2008J; 
iKhochfar fc Ostrikerll2008r ). The second mechanism, based 
on the inflow within the disc, is determined by the self- 
regulated disc_instability, and is the topic of a companion 
paper by ICacciato et al.l (|2011T ). The third, based on mass 
streaming in from outside the disc, represents an external 
source of energy that is determined by the cosmic growth of 
structure and is independent of the disc instability. In this 
work we focus on this external mechanism. 

So far, different authors seem to have reached dif- 
ferent conclusions concerning the possible role played 
by these 'gravitati o nal heating' mechanisms at z ~ 2 
(IGenzel et al.l 120081; IKhochfar fc sllkl 120091; iLehnert et al 



20QC ; lElmegreen fc Burkerj |2010| ; iKlessen fc Hennebelle 
20ld ; IKrumholz fc Burkerdl201Ch . This uncertainty natu- 



rally stems from the fact that the power provided by cos- 
mological in-streaming is in the same ball park as the tur- 
bulence dissipation rate. This motivates the more detailed 
analysis presented in this paper. 

The evolution of the gas mass in a disc galaxy can be 
described as a quasi-steady-state solution of a simple dif- 
feren t ial equation of mass conservation dFinlator fc Dave 
2008; Dek el et all l2009bl ; iBouche et al"1l20ld ; iDutton et al. 
l20ld ; iDave et alll2011al lbl. see below). In turn, the gener- 
ation of turbulence and its dissipation in a steady state 
is gover ned by an analogous equation of energy conserva- 
tion (e.g. IKhochfar fc Silk1l2009l ; lElmegreen fc BurkertfeOld ; 
IKlessen fc HennebelldbOlOT ). Together, these equations help 
constrain the disc instability, as the Q parameter depends 
on disc mass and turbulence. 

There is numerical and observational evidence for 
marginal instability with Q ~ 1 in disc galaxie s. Simulations 
reveal Q ~ 1 in high-redshift, gas-rich discs dlmmeli et al 



2004 ICeverino et~aT]|20ld ; iGenel et al]|2012l ; iHopkins et a! 



20111 ) and Q 



3 in low-redshift stellar-dominated 



tract i on or tneir energy ana momentum irJirnbon n a, ueken zuua; ivvestian et auizuin; ivan qer riruit <; 

2003 |; lKeres et allfflOarPekel fc Birnboiml200d ; lKeres et all lYim et alll201ll; IWatson et al.l I2012T ) and 



discs (|Hoh]||l97ll;lAthanassoula fc Sellwoodlll98d ; iBottemal 
120031 ; iMartig et all I2009T ) . Similar estimates are obtained 
from observed discs, both in local galaxies dLerov et all 
20081 ; IWestfall et"aTI l201ll; Ivan der Kruit fc Freeman! l201ll 



200d : iDekelet all l2009al l At low redshift, the flow con 
sists of c old gas clouds that 'rain' from the hot halo onto 
the disc JMaller fc Bullockl |2004 IDekel fc Birnboiml 120081 ; 
iKeres fc Hernauistll200d ) . Such flows may be efficient ly con- 
verted into turbulence (JKlessen fc Hennebellel |201y) , pro- 
vided that the in- streaming gas has comparable density to 
the disc ISM (e.g. IDekel etal]|2009bl ). 

The role played by stellar feedback in driving the ISM 
turbulence is controversial. On one hand, theoretical esti- 
mates and numerical simulations argue that stellar feed- 



z ~ 2 discs 
(|Genzel et al.ll2011f ). However, the parameter Q applies to 
linear perturbation theory and caution is required in com- 
paring it with estimates from the high non-linearly evolved 
(observed and simulated) galaxy discs. 

The instability of a disc and the associated level of 
turb ulence may be coup led via a self-regulation mechanism 
(e.g. jDekel et al.ll2009bl ). If the turbulence is driven by inter- 
nal processes that themselves depend on the disc instability, 
such as star formation or internal torques that cause mass 
inflow, the system may relax into a steady-state by a self- 
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regulation loop. In this case the disc maintains marginal 
instability, Q ~ 1, as the turbulence adjusts itself to the 
proper value dictated by the gas surfac e density and the an- 
gular v elocity. In a companion pap er (|Cacciato et al.ll201ll . 
see also iKrumholz fc BurkertluOlOT ). we impose Q = 1, and 
analyse the steady state solution of the mass and energy 
equations under the assumption that the energy source for 
driving the turbulence is the inflow down the potential gra- 
dient within the disc. The rate of this inflow adjusts itself to 
compensate for the turb ulence dissipative loss es such that 
Q = 1 is maintained. In ICacciato et al.l (|2011T ) we address 
the instability of a two-component disc, with gas and stars of 
different velocity dispersions that gradually exchange mass. 
We find there that discs tend to 'stabilize' at low redshift as 
the disc becomes d ominated by the 'hot' stellar component. 
iForbes et al.l (|2011f ) study a similar scenario including radial 
variations within the disc. 

In this paper, we study the steady-state solution of sim- 
ilar mass and energy conservation equations, but focus on 
an external energy source, carried by the cosmological in- 
streaming of gas. We address two main cases. First, a case 
where the system is governed by the external source alone 
and the instability is not self-regulated, and where the ef- 
ficiencies of the various physical processes are constant in 
time. Second, an alternative case where the instability is 
self-regulated by a duty-cycle for instability and star forma- 
tion. 

This paper is organized as follows. In Section [2] we 
present the equations of mass and energy conservation and 
their steady-state equations, and introduce our parameter- 
ization of the relevant physical scenarios. In Section [3] we 
investigate the non-self-regulated case I, with fixed efficien- 
cies of the physical processes, and predict the evolution of 
Q and a/V TO t- In Section [4] we study the self-regulated case 
II, where we impose Q — 1 and introduce a duty cycle for 
instability and star formation. In Section [5] we put our re- 
sults in the context of other results from the literature. In 
Section [6] we conclude and discuss our results. 



2 EQUATIONS FOR MASS AND ENERGY 
CONSERVATION 

In this section we present the basic equations for conserva- 
tion of mass and turbulent energy in a gaseous galactic disc. 
Source and drain terms for gas mass and gas turbulent en- 
ergy are identified. As a result, the gas mass and gas velocity 
dispersion that characterizes the turbulence are computed in 
a steady-state solution as a function of the incoming supply 
rate of cosmological gas and the parameters that character- 
ize the various relevant physical processes. 



and inflows inside the disc into the bulge. Eq. (|7|) has a 
very simple and instructive solution if the sink terms can 
be written as M s j n k = A/ g r _1 , i.e. if they are proportional 
to thegas mass itself with a 'sink timescale' proportionality 
factoiQ If Mcosmo and r vary on a timescale longer than r, 
the solution is 



M, 



-t/T 



-t/T 



(8) 
M g M COBmo e- I/T . (9) 

For t 2> t, it reduces to a steady state solution with 

Mink ~ Mcosmo (10) 

M g « (11) 

M g « M cosmo T, (12) 

in which the sink term M s i n k ad justs itself to matc h the 
external source term Mcosmo fsee lBouche et al.ll2010f ). The 
range of validity of this solution is discussed in detail in 
Appendix lAl 

Assuming that the disc has reached the steady state, we 
use the results from mass and energy conservation to derive 
the turbulent velocity as follows. We start by considering, 
for simplicity, only the star formation (SF) part in the gas 
mass sink term, in the form 



Ms 



Mg 

tSF 



and obtain from mass conservation 

MsF = Mcosmo , 

and 



M e = M C( 



3 £sf- 



(13) 



(14) 



(15) 



In analogy, for the turbulent energy E tur b we consider 
a sink term in the form of a dissipation rate, 



Edis = 



-Eturb 
tdis 



This leads in steady state, in analogy to eq. (|14|l . to 

-t-'dis — -LVcosmo j 



(16) 



(17) 



where £ cosmo is the rate of in-streaming energy. The ana- 
logue to eq. (|15|l is then 



-C'turb — -C'cosmot'dis 



(18) 



We approximate E tur b ~ M g a 2 and E coslno m M coamo V^, 
where the in-streaming velocity Vm, as well as the rotational 
velocity Kot, are a ssumed to be comp arable to the virial 
velocity of the halo (jDekel et al. 2009a), and the conversion 
of in-streaming kinetic energy to disc turbulence is assumed 
to be efficient. Eqs. ([15} and (fl8|) then yield 



2.1 The backbone steady state model 

The basic equation for the gas mass budget of a galactic disc 

is 



M E = M c , 



Msin 



ink i 



(7) 



where M g is the disc gas mass, M coamo is the external source 
term that represents the cosmological gas inflow rate, and 
M s i n k is the sum of different kinds of 'sinks' that empty the 
disc of its gas, including star formation, galactic outflows, 



(19) 

A very interesting feature of eq. ()19p is that a/V TO t turns 
out to be independent of the cosmological accretion rate it- 
self. This unique feature stems from the facts that (a) in 
steady state both the sinks of SFR and dissipation rate ad- 
just themselves to the corresponding supply rates, (b) in our 



1 As we discuss later, t is related to the dynamical time of the 
disc. 
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current model, the cosmological supply is a common source 
for both the mass and turbulent energy of the disc, such that 
the cosmological input always provides the same specific tur- 
bulent energy, and (c) this supply is determined externally, 
independent of the conditions in the disc. The implication of 
this special property of our current model is that a/V TO t is 
expected to be invariant under variations in the mass input 
rate, which is probably the main source of variation in the 
galaxy properties related to disc instability, both between 
different galaxies and as a function of time in the history of 
each individual galaxy. 

Both timescales in eq. ()19[) are expected to be related 
to the dynamical time of the disc, 



tdis — 'Ydist'dyn ■ 



(20) 



with 7di s a constant parameter with a likely va lue in the 
range 1-3 (|Mac Low et al.lll998l ; rGammiell200ll ). and 



tsF 



£dyn 
CSF 



(21) 



with e SF ~ 0.02 (|Silklll997l : rGenzel et al.ll2010l ). With these 
fiducial values eq. (|19|l becomes 



77— = V £ SF7dis « 0.2 . 
Vrot 



(22) 



In addition to the uncertainty in these parameters, several 
numerical factors of order unity have been omitted in this 
simple derivation, which will be recovered in Section 12.31 
Despite the fact that the numerical values of esf, 7dia and 
the other parameters are not known to great accuracy, we 
learn from eq. i)22[l that the available power in the external 
accretion is in the same ballpark as the power required at 
z ~ 2 for maintaining the turbulence in the discs. As long 
as the conversion efficiency of that energy into turbulent 
energy is high, this could in principle be the main driver of 
disc turbulence. At low redshift, the cosmological accretion 
carries more than enough energy to maintain the observed 
cr/V r ot « 0.05. 

As an aside, the contribution of stellar feedback to driv- 
ing turbulence in the disc can be estimated in a similar 
way, replacing the gravitational potential V in by the energy 
provided by stars per unit stellar mass formed, Vp B . Thus 
-Ecosmo ~ M cosmo V r 2 ot is replaced by Efb = Msf • Vp B - With 
~ 10 erg released by each supernova and one supernova 
per 100 Mq of stars formed, the released energy corresponds 
to Vfb ~ 700 kms -1 , well above what is needed for driv- 
ing the observed turbulence. However, the vast majority of 
the energy emitted is ra diated away in the initial phases of 
the s upernova evolution (jDekel fc Silklll986l ; lThornton et al.l 
1 19981 ). such that the energy available to be deposited in the 
ISM corresponds to only Vfb ~ 100 kms -1 , lower than the 
required energy. Furthermore, a large fraction of that energy 
is likely to drive out flows from the disc rather than turbu- 
lence inside the disc (|Mac Low fc Ferrarall 19991 ; Ijoung et al.l 
120091 ; lOstriker fc Shettvil201ll ). 



2.2 Mass steady state 

We now generalize the mass sink term in eq. (Q to include 
several additional processes, and have at steady state 



M s = M c , 



MsF — M w out + My, in — A/in 



0. 



(23) 



In the following, we describe the parameterization of the 
various sink terms. Table [l] gives an overview of all model 
parameters. 



The star-formation rate (SFR) is defined as 



M SF = De S F 



M^ 

tdyn 



(24) 



where D is a duty cycle (with a fiducial value of 1, to be 
discussed in Section I4.1|l and esf is the fraction of the gas 
that turns into stars every dynamical time idyn whenever 
star formation is 'on' (D = 1), as in eq. (12111 , 

• It is assumed that stellar feedback blows galactic winds 
at a rate that is proportional to the SFR with a mass-loading 
factor rj, namely 



M„ 



V M S f, 



(25) 



where typically r\ ~ 1. Galactic winds can also be 
dri ven by feedback from a n Active Galactic Nucleus (AGN; 
e.g. iNesvadba et al.ll2008l ). which, for simplicity, we include 
in the parameter 77, since AGN activity and star formation 
are often concurrent. 

• It is assumed that some fraction 7 rc i of the mass that 
was blown out into the wind comes back as an instantaneous 
fountain: 



M,, 



7rcl • 7/ ■ M S F ■ 



(26) 



The fiducial value we choose is 7 rc i = 0, namely no returning 
winds. 

• Disc instability is associated with torques that d rive 
angular momentum out and mass in (|Gammiel 120011 ) , to 
form a bulge at the disc centre. Part of this is clump migra- 
tion, whose rate could be computed by dynamical friction or 
clump-clump interactions. The inflow rate can be expressed 
a.s 



Mini = D7inf— - = 7inf£ S F lM SF ; 
Idyn 



(27) 



where 7j n f is the inflow 'efficiency' per dynamical time 
JDekelet al.l l2009bl . eqs. (19) and (24)). Based on eq. ©, 
we estimate 7i n f ~ 0.02 at z ~ 2, and 7i n f ~ 0.001 at z ~ 0. 
To avoid underestimating the effect of disc inflows, we use 
a fiducial value of 7infeg" F = 1. The parameter D in eq. (|27|l 
is the same duty cycle as in eq. (|24|l for the SFR, assumed 
to be determined by the instability duty cycle. 



M e 



By solving eq. (|23|) . we obtain the steady-state solution 



(28) 



-'Wcosmotdyn 



£>6Sf[1 + »?(1 — 7rcl) + 7inf£sF] 



The contribution of each term is easy to understand qual- 
itatively. The gas mass is proportional to the cosmological 
supply rate of gas. More vigorous outflows (large rf) reduce 
the gas mass in the disc, unless they are largely recycled 
(large 7 rc i)- A duty cycle with longer 'off' phases (small D) 
leaves more gas in the disc. A higher efficiency of star for- 
mation (large esf) reduces the gas mass. A stronger inflow 
inside the disc (large 7i n f) also reduces the disc gas mass. 
Using eq. (|24|l . we also obtain 



Msf = 



M c , 



1 + r?(l - 7 rc i) + 7infe s 



(29) 



Table 1. Model parameters. 
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Parameter Fiducial Value Definition 



CSF 


0.02 


D 


1 


>l 


1 


7rcl 





7inf 


0.02 


7dis 


2 


(i 


1 


vm 





U 


I/a/2 


J- 1 






Equation 



Star-formation efficiency per disc dynamical time i ]24[ ) 
Duty cycle for instability, star- formation and inflows H24IU27H 

Wind mass-loading factor 1251 

Fraction of instantaneously recycled wind H26II 

Fraction of gas inflowing inside the disc per dynamical time H27H 

Ratio of the turbulence dissipation timescale to the dynamical time (1201 ) 

Fraction of the in-streaming kinetic energy that turns into turbulence J30I ) 

Fraction of the disc inflow potential energy that turns into turbulence J30I ) 

Ratio of V ro t to the in-streaming velocity 1311 1 

The slope of the logarithmic rotation curve, (log Vrot/logr) <£3j 



2.3 Turbulent Energy steady state 

Gravitational energy can be transferred into turbulent en- 
ergy to compensate for the dissipative losses in two gen- 
eral ways. If a large fraction of the incoming streams 
that hit the disc is in dense gas clumps, they can trans- 
fer momentum into the disc gas, an d thus convert the 
stream kinetic ene r gy into turbulence IIGenzel et al.l 120081 ; 
iDekel et al.1 l2009bl ; lElmegreen fc BurkertJ 120101) . Alterna- 



tively, the rotational and potential energy of the disc mass 
becomes available due to the instability-driven mass inflow 
within the disc, including clump migration, where turbu- 
lence is generated by the same torques that are respon- 
sible for angular-momentum outflow and the associated 



mass inflow JDekel et al.ll2009bl; iKrumholz fc Burkertj[201ol : 



ICacciato et al.ll201ll ; iForbes et al.ll201lf ). The quasi-steady 
state of the turbulent energy in the disc is then described 
by 



K 



9 = Si(M 



+ Mw,m)0.5^ + UMnf Kot 



(30) 



(Msf + M w ,out + M inf )1.5a 2 



M g 1.5a 2 t^ s = 0, 



We solve eq. $U} for a/V T0 t, using eqs. flUJ and (J2H1), 
and obtain 

<? 2 = 2/3 fr (1 + q + TinfggF 1 ) + g m 7inf 6^(^/0.5) 
V? ot k 2 /0.5 (esFTdis^)" 1 + (l + i) + 7inf£sp) 

In the limftl 2 ] 1 + rj + 7inf£gp <C (esF7dis-D) _1 and with u 2 — 
0.5, we obtain 



77 — ~ 0.16< v /esF,o.027dia,2-D 

Vrot 



(33) 



X \J &(l + T) + ^ni^SF 1 ) + Cm^nf^SF 1 > 



where esF,o.02 = esf/0.02 and 7 d i s ,2 = 7dis/2. 

As already noted for the approximate steady-state solu- 
tion in eq. Q22p . the more detailed result for a/V ro t in eq. (|32|) 
is also independent of M cosmo , because in our current model 
the latter controls both the incoming energy and the disc gas 
mass that is involved in the turbulence, so the varying cos- 
mological supply always provides the same specific turbulent 
energy. 



where the different terms are as follows: 

• It is assumed that both the cosmological accretion and 
the recycled wind arrive to the disc at a speed 



Vin = Kot/w, 



(31) 



and that a fraction £; of the kinetic energy they carry is 
converted into disc turbulent energy. We take u 2 — 0.5 (see 
Appendix[Cj, and a fiducial maximum value fj = 1, so that 
we can examine the maximum possible contribution of in- 
streaming energy. 

• It is assumed that as gas flows in or migrates to the 
galaxy centre, potential energy, which is released at a rate 
of ~ Mj n fV^o t , is transformed into turbulent energy with 
an efficiency £ m . See Appendix [B] for a discussion of this 
assumption. We later examine both the effects of £ m = 
and £ m = 1. 

• The gas sink terms, i.e. star formation, outflowing winds 
and inflows to the bulge, take their share of turbulent energy 
when they leave the disc. 

• The turbulent energy dissipates on a dissipation 
timescale idis = 7di s idyn, as in eq. (|20[) , with a fiducial value 
7dis = 2. 



3 CASE I: NON-REGULATED DISC 
INSTABILITY 

3.1 The relative roles of in-streaming and disc 
inflow 

The approximate expression -Ecosmo ~ M cosmo V^o t used for 
the simple derivation in Section \2. II does not distinguish be- 
tween cosmological in-streaming kinetic energy and gravi- 
tational potential energy released during inflows inside the 
disc. Each of these components carries similar specific energy 
of ~ V^ot (Appendix [C| . Inspection of eq. I)33p elucidates 
their individual roles. 

The conversion of rotational and potential energy into 
turbulent energy depends on inflow in the disc, i.e. £ m cou- 
ples to 7i n f . The ability of the kinetic energy of in-streaming 
mass to be converted into turbulence in the disc depends on 
complex physical processes in the vicinity of the disc, which 



With the fiducial values (Table [TJ, 1 + r\ 4- 7i n f£g 



3 and 



( e SF7dis.D) = 25, so we consider the approximation 1 + r\ + 
7inf e sp *C (eSF7dis£ ) ) -1 a good one. More generally, 77 and 
7inf e gp are n °t expected to have a value of more than a few, 
while (esF7dis£ ) ) -1 ' s expected to have a value of several tens. 
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we simply parameterize with £j. If & > 0, the accretion 
energy acts as a direct external driver of turbulence inde- 
pendently of disc inflows or winds. However, the presence of 
disc inflows and/or winds enhances the contribution of in- 
streaming in driving turbulence, i.e. £j couples both to 7i n f 
and rj. Thus, the relative contribution of kinetic and poten- 
tial energy of cosmological origin depends not only on the 
intrinsic efficiencies £; and £ m , but also on the importance 
of winds and disc inflows. 



3.3 Evolution of Q 

If the values of the parameters D, £j, and £ m are fixed, both 
M s and a are determined by a balance between the exter- 
nally set in-streaming, the star formation, the winds, and 
the inflows in the disc. In this case, the Toomre Q parame- 
ter is not necessarily locked to Q ~ 1. Self-regulation of the 
instability at Q ~ 1 requires that the relevant physical pro- 
cesses, such as the inflow rate in the disc, the SFR and the 
outflow rate, adjust themselves to maintain Q ~ 1, and this 
case, where the relevant parameters are not fixed, is deferred 
to Section [4] 

The Toomre stability parameter is 



3.2 The effects of galactic winds and disc inflow 

By comparing eqs. (|28[) and (|15[) we can verify the role played 
by winds and inflows in the disc on the gas mass and hence 
the SFR. The gas mass and SFR are suppressed by escaping 
winds. In the limit of very strong winds and little recycling, 
77(1 — 7 rc i) 2> 1 + 7i n f £g F , the gas mass and SFR are roughly 
proportional to rj. If the winds are fully recycled, 7 rc i = 1, 
then they have no net effect on the gas mass and SFR, given 
that recycling is assumed to be instantaneous. 

On the other hand, the value of a is independent of the 
recycling rate 7 rc i. This is despite the fact that the origin 
of the dependence on rj does depend on j TC i. For example, 
If there is no recycling, 7 rc i ~ 0, a is higher for larger r) 
because the gas mass is smaller while the energy input is 
the same. On the other hand, if the outflows are all recycled 
back to the galaxy, 7 rc i ~ 1, a is higher for larger rj because 
the returning winds add to the energy provided by the cos- 
mological accretion while the gas mass remains unchanged. 
We note that, regardless of the level of recycling, stronger 
outflows drive an increase in a, despite the fact that we have 
ignored direct local deposit of feedback energy in the disc 
gas. This is done indirectly, either by lowering the gas mass 
or by adding to the energy brought by external accretion. 

The inflow in the disc depletes the gas mass and thus 
suppresses the SFR, via 7i n f£sF' a quantity of order unity. 
The turbulent velocity a is even more sensitive to the in- 
flow in the disc because (a) when the gas mass is depleted 
the same amount of energy input by in-streaming results 
in higher a (via £i), and (b) the inflow down the potential 
gradient in the disc contributes energy to driving a up (via 

Cm). 

Can the observed velocity dispersion at z ~ 2 be pri- 
marily driven by the cosmological in-streaming? Recall that 
the observed values are typically cr/Kot ~ 0.2, perhaps even 
~ 0.3. If the conversion efficiency of in-streaming energy 
to turbulence is high, fj » 1, and the SFR and dissipa- 
tion rate are at their fiducial values, the term referring 
to the streams by themselves already provides cr/Kot ~ 
0.16, which is in the ballpark of the desired value, though 
slightly short. With the fiducial depletion by winds, r) ~ 1, 
this becomes cr/Kot ~ 0.2. With the fiducial disc inflow, 
7inf£gF ~ 1 (and £ m ~ 1), even without winds, it becomes 
cr/Kot ~ 0.3. Adding the fiducial winds and disc inflow, we 
obtain cr/Kot ~ 0.3. We conclude that if somehow £, ~ 1, 
the gravitational energy associated with the inflow, and in 
particular the clumpy gas streaming into the disc, can have 
a significant contribution to the disc turbulence, which can 
be naturally aided by outflow depletion and disc inflow. 



Q = 



7rGE E 



(34) 



where n = y/2vtt. Using fi = V rot /R, Kot = GM tot /R, 
Mtot = Mg/6, and M g = -k R E g , we obtain 



Vrot 



(35) 



Then using the approximate solution in eq. ((33} for cr/Kot, 
and assuming a flat rotation curve v — 1, we obtain 



Q ss 0.68<5 33 \/esF,o.027dis,2-D 

X y£i(l + J7 + 7mf£sp) +£m7infesF> 



(36) 



where the disc mass fraction (eq. (JSJ) is expressed by <5o.33 = 
5/0.33. Substituting the fiducial values from Table [TJ we 
obtain Q w 1.18<5^~3 3 . 

From eq. (|35[) . in the non-self-regulated case studied 
here, where cr/Kot is constant, we learn that Q scales with 
mass and time as Q oc S^ 1 , where 6 — M g /M to t within the 
disc radius. We can evaluate the time evolution of Q through 
5 as follows. From eq. (|28[) . the steady-state solution is M s oc 
Afcosmoidyn- Assuming that the disc radius is proportional 
to the halo virial radius, R = A_R v ir, with A a constant spin 
parameteqj, we obtain idyn oc iHubbio, which in the Einstein- 
deSitter phase that is approximately valid at 2 > 1 gives 



tdynOc(l + z) 3/2 . 



(37) 



Based on the EPS approximation (confirmed by fits to cos- 
mological simulations), the cosmological input rate can be 
approximated in the Einstein-deSitter phase by 



Mc, 



OC Mvir (1 + z) 



5/2 



(38) 



(|Neistein et all I2006J ; iNeistein fc Dekell 120081 ). Assuming 
that the halo mass profile is roughly M(r) oc r, we have 



3 This is based on the assumption of gas angular moment um con- 
servation during the collapse to the disk IIMo et alj|l998» . While 
numerical st udies have suggeste d that this picture neglects many 
details (e.g. IZavala et alj I2008J; button fc van den Boschl [20091; 



iDanovich et al.ll2012l ; lAumer fc White! |20~ V it provides a fair 



matc h to observed disk size s (e.g 
200l|: iMaller fc Dekell |2002|; lDut„ 
2010t button fc van den Boschll2012l) 



Mo et al. 
tton et alj 120071 : 



iBullock et al 



iBurkert et al 



Cosmological inflow and disc turbulence 7 



Mtot — AM v i,Q The above finally yield for the non-self- 
regulated model 

QocfT 1 oc (1 + zy 1 . (39) 

We learn from eqs. (|36p and ()39p that the discs tend to be 
unstable at high redshifts, Q < 1, and then evolve toward 
stabilisation at later times, Q > 1. We can see that the 
growth of Q in time is driven by the decline of the cosmolog- 
ical accretion rate, eq. (|38|l , being steeper than the increase 
in time of idyn- This trend is in at least qualitative agreement 
with observation s that find Q to be of order a few in local 
disc galaxies (e.g. lLerov et al.12 008; Wa tson et al.l2012T). but 
aroun d, or even below, unity at high redshift 1 Genzel et al.l 

[201I]). 

Eq. ()36[) shows that at high redshift, violent disc insta- 
bility with Q < 1 is naturally driven by the high 5, which 
results from the intense in-streaming rate, eq. ()28|) . The in- 
streaming power, even when the efficiency for driving tur- 
bulence is high, £j = 1, may not be able to drive turbulence 
with cr/Kot high enough for balancing the high S and thus 
stabilising the disc. The instability is associated with disc 
inflow of a large 7i n f, which can drive further turbulence 
both by providing energy that is converted to turbulence 
if Cm is high and by depleting the disc mass (see Section 
[4.30 . High-redshift ob servations reveal un stable discs where 
Q ,$ 1, with a high 3 IITacconi et al.ll201C) as w ell as a high 
r\ (e.g. ISteidel et alll2010l : iGenzel et aljboilf ). The latter 
may indicate that feedback also plays a role in maintain- 
ing Q ~ 1, either indirectly by depletion of the disc gas 
through outflows, or by direct injection of feedback energy 
into turbulence (not considered here). Some combination of 
disc inflows and feedback/outflows is capable of providing 
turbulence at the level that would stabilize the disc, and 
this may occur in a self-regulated way that keeps the disc 
marginally unstable. At low redshift, where 5 <C 1, Q could 
be higher than unity even without the presence of winds 
and disc inflows (see also our complementary calculation in 
ICacciato et al.ll201ll 'l. 

We can replace the assumption of a constant loading 
facto r rf with a prediction for momentum-driven winds, r\ oc 
U _1 (jMurrav et al.ll2005f h Using the standard cosmological 
virial relation this implies 

(40) 



Qocikr 1/2 (l + z) 



-7/4 



(42) 



r,ocM~ 1/3 (l + z)- 1/2 . 



Plugging this in eq. (|33|l . in the limit of strong outflows 
r\ 3> 1, we obtain 



-1/2 



v n 



oc7 ? 1/2 ocI/~ i/ «M- i/0 (f + z) 



1/6/ 



-1/4 



(41) 



Then from eq. (|35[) . using eq. (|28[1 that implies S oc rj 1 (1 + 
z), we obtairM 



4 This is formally exact for an isothermal sphere, which is a rea- 
sonable ap proximation for the dark matter halo. Since typically 
A ~ 0.03 IjBullock et al.l 120011 ). and the baryons in the central 
galaxy are known to be les s massive than fS 0.03M v ; r I IGuo et alj 
2010; Bcfir oozi et alj|2010f ). we find that this is a reasonable ap- 
proximation altogether. 

5 Note that there could be tension between the limit 1 + w + 
7inf e sp ^C (esF7dis£ ) ) -1 hi which eq. 13311 is valid, and the limit 
n 2> 1 in which eq. fl41| l is valid. They may both be valid if rj ~ 5, 
say. If n is larger, we can relax the first condition and use eq. 1321 1 



We learn that the variation of n according to momentum- 
driven winds enhances the redshift dependence of Q com- 
pared to the constant r\ case, eq. (|39p . Unfortunately, the 
slight growth of a/V I0 t in time in eq. (|41|) makes the agree- 
ment with the observed decline of a/Kot even worse than 
in the constant n case. The consideration of momentum- 
driven winds also introduces a mass dependence. The pre- 
dicted weak mass dependence of cr/Kot oc (1 + z)~ 1 ' 4 
seems reasonable, although th e observational trend is not yet 
well established. For example, iKlessen fc Hennebellej (J2010h 
find that ff/Kot oc M~ x ' 3 in the local Universe, while 
Ivan der Kruit fc Freeman! (|201ll ) report a/V TO t «const. At 
z ~ 2, the largest existing observational sample shows no 
clear trend, though still wi th significant scatte r that makes 
the situation inconclusive (|Mancini et al.ll2011T ). 

The observed mass dependen ce of Q in galax y discs 
is not yet well established. While iDalcanton et alj (J2004 ) 
find a shar p threshold for the on set of instability at Vrot > 
120 kms"\ IWatson et ail (|2012T l find no trend of Q with 
galaxy mass. Our predicted mass dependence of Q depends 
on the scaling of 8 and the different model parameters with 
mass. In the case of constant model parameters, we obtain 
that Q is independent of mass, eq. (|39[) . and in the case of 
momentum-driven winds scaling, we obtain that more mas- 
sive discs are less stable, eq. (1421 1. 



4 CASE II: SELF-REGULATED DISC 
INSTABILITY 

In Section [3] we have shown that case I produces gas discs 
in which cr/K-ot is constant in time and Q is gradually in- 
creasing to values larger than unity. While the systematic 
increase in Q toward low redshifts may be consistent with 
the observational trend toward a larger abundance of stable 
discs, the constancy of a/V ro t is clearly in conflict with the 
observed decline of this quantity. In this section, we appeal 
to an alternative case II, where self-regulation at Q ~ 1 is 
imposed at all times. We address the possibility that the self- 
regulation is achieved by periodic episodes where the insta- 
bility and the associated SFR and disc inflows are 'on' or 'off 
with a duty cycle D < 1 (see also iMartig et al.ll2009T ). Self- 
regulation may alternatively be achieved by adjustments of 
the dis c inflow rate, via 7i n f and £ m , as in our companion 
paper (|Cacciato et al.ll201lr ). We also check the effect of a 
systematic decline with time of the accretion conversion effi- 
ciency £i , perhaps reflecting the evolution from dense narrow 
streams to a wide-angle accretion and a gradual decrease in 
the gas fraction and degree of dumpiness in the accreting 
gas. 



instead of eq. 1331 1. We then find that the scaling of Q in the 
r\ ^> 1 limit is bounded by the scalings in eqs. 1391 1 and 1421 1. 
and the scaling of cr/Vrot is similarly bounded by the flat value 
of eq. 1321 and the weak variation in eq. 1411 1 . 
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4.1 Maximum conversion efficiencies and a duty 
cycle 

We assume here that once Q rises to slightly above unity, the 
disc tends toward stabilisation, star formati on is suppressed 
(|Kennicuttl 1989; M artin fc Kennicuttll200ll ). and the galac- 
tic outflow as well as inflow within the disc become weaker 
too. This allows E g to increase in response to the continuing 
cosmological accretion, while the growth of a slows down. As 
a result, Q tends to decrease to slightly below unity, back 
to an unstable phase where star formation, outflows and 
disc inflow resume, and so on. We thus expect oscillations 
about Q ~ 1. We model this by a duty cycle D < 1 that 
represents the fraction of time when the instability is 'on'. 
Here we adopt maximum efficiencies for energy conversion 
to turbulence, £i = £ m = 1. From eq. (|36|) . with Q = 1, we 
obtain 



D 



2.2<5q.33 



esF,o.027dis,2(l + V + 27 inf e SP ) 



(43) 



With the fiducial parameter values (Table[T} , and a high 
cold mass fraction characteristic of high-redshift discs, 8 ~ 
0.33, we have D £, 1. When the gas fraction is even higher, 
the outflow rate is lower, or the disc inflow rate is lower, 
eq. (|43[) may give D > 1, which is clearly unphysical. This 
is another representation of the result from Section [33] that 
at high redshift, when 8 is high, in the absence of winds and 
when the disc inflow is ignored, there is hardly enough power 
in the cosmological in-streaming by itself to drive the high 
turbulence required for Q — 1. The fiducial winds, r\ ~ 1, 
or the fiducial disc inflow, 7i n feg F ~ 1, help obtaining a 
physical result with D < 1. Another way to obtain D < 1 
may be if during the violent instability phase the timescales 
for star formation and inflow somehow become shorter than 
the disc dynamical time idyn, while the dissipation timescale 
is long compared to idyn- 

On the other hand, at low redshift, with 8 < 0.1, and 
even more so if 7i n f and r\ are non-negligible, the resulting 
duty cycle D drops significantly below unity to ,6 0.1 in 
order to keep Q = 1. From eq. (|43[1 . the scaling of D with 
mass and redshift is Doc M°(l + z) 2 ' 3 in the limiting case 
of weak winds (r) < 1), and D oc M 1/3 (1 + z) 7/6 in the limit 
of strong momentum-driven winds with no recycling (rj 3> 1 
and 7 r d = 0). fn the general case between these limits, D is 
declining with time close to linearly with (1 + z), such that 
at later times one expects to detect a smaller fraction of the 
galaxies in the unstable phase. 

When forcing Q — 1 and keeping D free, from eqs. (|28[) 
and (|33|l and the definitions for k and E g , we obtain 



7disM COBmo (l + rj + 27 inf e S p) 

7T y 3y/2u 2 G 2 (l + Tj(l - 7rcl) + TM^Sf) 

1 
_ 2 /35Myr\ ( M cosmo 



(44) 



and 



36M Q pc" 



GM 



1 M Q yr- 



3 7dis(l + V + 27 inf e S p) 



3V2u 2 (l + r)(l - 7rcl ) + 7infesF) 

1 

M 

■"■^cosmo 



(45) 



where the second equality in each of these equations is calcu- 
lated using our fiducial parameters from Table [l] In contrast 
with the solution in Section [3j here a/V TO t does depend on 
the external cosmological in-streaming rate M cosmo . This is 
a result of the imposed self-regulation, where a decrease in 
Af cosmo induces similar decreases in a /V YO % and in 8 such that 
Q remains constant, eq. (|35[l . These evolution trends agree 
with the observed trends better than those predicted in case 
I. In particular, the results obtained for M cosmo = 1 Mq yr _1 
and k = (35Myr) , as roughly appropriate for the Milky 
Way, are close to the typical values in local disc galaxies, 
and if M cosmo is scaled up to « 100 Mq yr _1 , as appropriate 
for high redshift discs, one obtains E g w 200 Mq pc~ 2 and 
a ~ 80kms~ , which are indeed in the observe d ballpark 
(JForster Schreiber et af1l2009l : lGenzel et all201Cl ). In case II 
the gas mass M g declines at a slower rate, because the overall 
gas depletion by star formation, outflows and disc inflow is 
slower in accord with the low duty cycle D. An examination 
of eq. (|28[) using the aforementioned scaling of D in this so- 
lution, as well as eqs. (|37[1 and (|38[1 . gives M g oc M(l-l-z) 1 ' 3 . 
Similarly, eqs. {44]) and (|45j) give E g oc M 1/3 (l + z) 7/3 and 
a oc M 1/3 (l + z) 5/6 , respectively, and a/V m t oc Af°(l + z) 1/3 , 
for any constant rj. These trends with redshift agree reason- 
ably well with observations as well. The associated evolution 
rate of the gas fraction, M s /M oc (1+z) 1 ' 3 , is however some- 
what weaker than observations in dicate for mass i ve ga laxies 
between z w 2 and z = (see iBouche et al.l (J2010J) and 
references therein). 

We note, however, that a decreasing star formation duty 
cycle, in particular one that decreases as strongly as 3 2 , 
should manifest itself in observations. If the 'off' periods 
are long compared to the galactic dynamical time, D< 1 
will result in a large fraction (~ 1 — D) of disc galaxies hav- 
ing significant gas discs with no star-formation. There are 
indeed several examples of observed gas discs that have low 
levels of star-formation, po s sibly as a result of their stabil- 
ity (e.g. iMartig et all 120091 : iMacLachlan et alj|201lh . How- 
ever, for typical d isc ga laxies this is not the case. Indeed, 
iLerov et all (J2008I ) and IWatson et all IJ2012T ) find no clear 
correlation between star-formation efficiency and Q. On the 
other hand, if the 'off' periods are short compared to the 
galactic dynamical time, they will result in discs being in 
a mixed state of 'on' and 'off' in different regions on the 
disc. In a statistical sense, this will show in observations as 
a decrease of the mean star-formation efficiency with cos- 
mic time. Such a trend may be observed, but probably not 
to the extent that is predicted by our model (JGenzel et all 
l20ld : iDaddi et alllioiol : iKrumholz et al.ll2012h . Disc galax- 
ies are observed to have, on average, similar star-formation 
efficiencies per dynamical time at low and high redshift. This 
remains an open problem with regards to this scenario. 



17.5 kms 



lM Q yr- 



4.2 Sub-maximal efficiencies and no duty cycle 

We now allow only a fraction £i of the incoming kinetic 
energy to be converted into turbulent energy, and solve for 
the value of £i that is required to keep Q — 1. We start with 
the accretion energy as the only direct driver of turbulence, 
namely £ m = 0, and assume D — 1. From eq. (|36[) . with 
Q = 1, we obtain 
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6 



£SF,0.027dis,2(l + T) + 7inf£ SF ] 



(46) 



This is a very similar expression to the one found for D in 
eq. (|43|l . However, the physical situation is different in the 
two cases. We note that <5 evolves differently because in the 
D/l case M s is larger than in the current case, as shown 
in Section 14.1 I In the current solution, E g is the same as in 
case I, i.e. it is decreasing in proportion to A/ cosmo . but due 
to the requirement Q — 1, here a is forced to decrease as 
well. This is achieved by a decreasing conversion efficiency 
£; from in-streaming to turbulence. 

The alternative, more physically motivated assumption 
of efficient conversion of the energy associated with the disc 
inflow, £ m = I (Appendix [B}, leads to 



ii 



2.2^.337^,2 - 50 7 inf 
CSF,0.02(1 + ?7 + 7mf e ; 



(47) 



SF - 



In both of these £j < 1 cases, the scalings of M s , 5 and 
g with mass and time are the same as in case I. From Q = 1 
2^, which implies 



we obtain cr/Kot = 8/ 
(T/Vrot ocM°(l + Z) 1 
without winds, and 
cr/Kot ocAH(i + 2 ) L5 



(48) 



(49) 



with escaping strong winds. The evolution of cr/Kot in the 
cases where £j < 1, with or without winds, agrees better 
with observations. 

At lower redshift, as the cold disc fraction 5 declines, 
the condition Q = 1 with D = 1 forces the conversion effi- 
ciency of the in-streaming to turbulence £; to decline as well, 
£i oc 8 2 , as in eqs. (|46l) and (|47|) . This means £, oc (I + z) 2 
without winds and £j oc M(l + z) 3 ' 5 with escaping strong 
winds. There could be several reasons for such a decline of £j 
in time. First, the conversion factor could be a growing func- 
tion of the gas fraction in the disc, namely of 5. Second, it 
is likely to be growing with the gas fraction in the accreting 
baryons, which is declining with time. Third, the evolution 
of the input pattern from narrow, dense streams at high 
redshift to a wide-angle accretion at late times makes the 
penetr ation of cold gas into th e disc less efficient at late r 
times (|Dekel fc Birnboiml [20061 ; Ivan de Voort et all I20TH ). 
The lower gas density contrast between the streams and 
the disc is another reason for a smaller fj, especially if it 
is somehow associated with a lower degree of dumpiness 
in the streams, as the coupling between in-streaming and 
disc arises fr om in-streaming gas of comparable density to 
the disc gas (jDekel et al.ll2009bl ). However, except for the 
first reason that relates to the gas fraction in the disc, the 
value of f j is determined externally, independently of the in- 
stability state of the disc or its other properties. Therefore 
a self-regulation loop is not expected to drive the required 
£i oc 5 2 for Q ~ 1. In this case, it seems that the evolution 
of f j can only match the requirements for Q ~ I by some 
coincidence. 

4.3 Self-regulated inflow in the disc 

Assuming fixed values for £j, £ m and D, the inflow rate 
within the disc may adjust itself to compensate for the 
dissipative losses and maintain the instability at Q ~ 1. 



This is the basis of t he an alysis in iKrumholz 
20101) 



, Forbes et al 
ICacciato et al.1 (|201 



et al.l 

ion]). 



Burkertl 



|201lf ). and our companion paper 
With our fiducial choice of parame- 
ters, including outflows, this can be achieved quite naturally. 
In the absence of outflows, the disc inflow rate should be 
comparable to or somewhat higher than the SFR, 7^ > esF, 
depending on the exact valu es of the parameters &, £ m , 7di 



and esF (see a discussion of IKrumholz fc Burkertl (|20101 ) in 
Section[S]). We note that an enhanced instability that boosts 
up the disc inflow rate would also enhance the SFR and out- 
flow rate, which would help the self-regulation. However, if 
such a case is also accompanied by faster dissipation, then 
increased inflows may not be able to solve the problem. 
The details of the inflow within the disc, including gas and 
stars, clumps and off-clump material, is being investigated 
via cosmological hydrodynamical simulations (Cacciato et 
al., in preparation). This has additional important conse- 
quences regarding the issues of bulge gro wth and feeding 
central black holes (jBournaud et al.ll2011al lbl). 



5 DISCUSSION: COMPARISON WITH THE 
LITERATURE 

Several earlier studies evaluated the possible role of in- 
streaming and disc inflows in driving the observed ve- 
locity dispersion in z ~ 2 discs, reaching seemingly 
conflicting conclusions. The first estimates of this kind 
(|Forster Schreiber et afll2006l : lGenzel et alj|2008l ) suggested 
that it is plausible that there is enough energy in the cosmo- 
logical in-streaming, depending on the exact n umerical val- 
ues of s everal parameters in their equations. In lDekel et al.l 
1)200913 ) the in-streaming energy explicitly depends on the 
unknown small-scale dumpi ness of the stre a ms, a nd could 
therefore go either way. In ICacciato et al.l (|201ll ) we ob- 
tain cr/Kot ~ 0.2 whe n choosing favourable values for 
the relevant parameters. iKlessen fc Hennebelld (|2010T ) es- 
timated that there is enough energy in t he cosmological 
accretion, and IKrumholz fc Burkertl (|2010j) reached a sim- 
ilar conclusion when examining the gravitational poten- 
tial e nergy that is rele a sed d uring mass inflows inside the 
discs. iKhochfar fc Silk! (|2009T ) concluded that only 18% of 
the in-streaming energy is required to reproduce the ob - 
serv ed values. In apparen t cont rast, iLehnert et al.l (| 20091 ) 
and lElmegreen fc Burkertl (|20101 ) conclude that there is not 
enough in-streaming energy to account for z ~ 2 turbulence. 

As we show in Section [2] the value of cr/Kot depends 
basically on two timescales, one associated with the decay 
of turbulence (parameterized with 7dia), and the other asso- 
ciated with the gas mass conservation in the disc (related to 
esF, r\ and 7i n f). All of the aforementioned studies involved 
a turbulent energy balance that is very similar to the one we 
consider in this work. However, they differ from one another 
in the choices of the turbulent dissipation timescale. More- 
over, some of these models do not explicitly include a mass 
steady state condition, and differ in the implicit assumptions 
they make regarding the timescale associated with the mass 
equation. A closer inspection of the different assumptions 
made, which we perform next, reveals how they lead to the 
app arently conflic t ing co nclusions. 

iGenzel et al.1 (|2008l ) do not explicitly write a steady 
state equation for the gas mass conservation, but their treat- 
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ment is equivalent to the simple equations in Section[2l only 
with tsF crudely approximated by the specific rate of cos- 
mological accretion t acc (from the growth of dark matter 
haloes), hence not addressing the possible difference between 
the rates of star formation and accretion. They remain ag- 
nostic as for whether there is enough in-streaming energy for 
driving the observed a/V ro t ~ 0.2 — 0.3, due to uncertainties 
in numerical values of order unity. Nevertheless, their fidu- 
cial valueqj indicate a small value of a/V TO t ~ 0.07, imply- 
ing that an additional source of energy is required. Apart 
from the choice of parameter values, there is a significant 
difference between our models. As t acc declines wit h time 
faster than the dynamical time, iGenzel et al.l (|2OO80 natu- 
rally obtain that a/V ro t decreases with time (see eq. (|19|) '), 
In our model, M s follows M cosmo (eq. (J2HJ)) and becomes 
smaller with time, so that the effective timescale for specific 
ma ss and energy gain is constant (and related to idyn). With 
the lGenzel et al.l (|2008l ) implicit assumption £sf = t acc , the 
gas mass remains constant with redshift, since the star- 
formation timescale becomes longer at the same rate that 
the accretion diminishes. 

iKhochfar fc Silkl 1)20091 ) find that a conversion efficiency 
of only £i = 18% is required to obtain high a/V I0 t as ob- 
served. Their semi-analytical model does not explicitly in- 
clude a steady-state solution to the gas disc mass as we 
do here, but gas velocity dispersions are obtained in their 
model as a result of similar physical considerations to the 
ones made here in Section 12.11 Therefore, it is straight- 
forward to compare the two analyses and find that they 
require a lower conversion efficiency £j as a result of the 
higher star-formation efficiency esf (a in their notation) 
they assume for their porosity-driven star-formation model, 

esf « 0.02(cr/10 kms" 1 ) « 0.15 at z ~ 2. 

The set of equations developed in Elmegreen fc Burkertl 

1)20101 ) is very similar to our equations here. They do not con- 
sider the winds (77) and disc inflow (7inf) terms, while they 
do consider a constant conversion efficiency £j (e in their no- 
tation) from in-streaming kinetic ener gy to turbulent energy 
and n o star formation when Q > 1. Elmegreen fc Burkertl 
(J2010j) focus on the transient state before the quasi-steady 
state, an issue we do not address here. Nevertheless, they 
obtain that once the quasi-steady state is reached, the accre- 
tion cannot be the dominant source of turbulence, because it 
cannot provide enough turbulence to keep the disc at Q ~ 1. 
When we neglect winds and disc inflows, we arrive at a sim- 
ilar conclusion, unless we choose favourable parameter val- 
ues such as £sf = 0.03 and 7di s = 3, for which we obtain 
cr/Kot = 0.24 (see eq. (13311). An add itional difference be- 
tween the lElmegreen fc Burkertl (|2010T ) work and ours is that 
they take the dissipation timescale to be the perpendicular 
crossing time, such that 7dis = t pclp u} = (H/a)(TTGY> g /a) — 
(cr 2 /cr7rGEdiac)(7rGSg/(T) = 5, while we assume 7di s to be a 
constant ~ 1 — 3. As a result, their solution at low accre- 
tion rates (and low S) is not as over-stabilized as ours (but 
still more stable than a case with a higher accretion rate). 
However, they do not directly discuss the case where the ac- 
cretion is low enough that the disc is over-stabilized at late 



times, where we have shown that a decreasing £; or D may 
allo w for a late marginally - unstab le configuration. 

iKlessen fe Hennebelkj (|2010r ) conclude that there is 
more than enough accretion energy to drive the observed 
turbulence at high redshift. They make a simple compari- 
son between the energy input rate 0.5Mi n V-^ and the dissi- 
pation rate, and apply relevant numerical values in different 
situations to see which one is larger. In the case of clumpy 
discs at high redshift, it seems that they used unrealistic 
numerical values. Using appropriate values in their eq. (8), 
either for individual clumps or for the whole disc, yields 
that a high conversion efficiency is required, £j ss 1, rather 
than £s <C 1 according to their estimate. This would imply 
that there is barely enough energy in the in-streaming for 
being the sole driver of the turbulence, in agreement with 
our results. In particular, for the case of single clumps, the 
result in their eq. (16) should be compared with the SFR in 
a single 10 8 Mq clump, i.e. ~ IMqjt' 1 , not with that of 
the whole disc, 10 — 50 Mq yr" 1 . Such a proper comparison 
would give a conversion efficiency of ss 1. For the case of the 
whole disc, their choice of 03D = 30kms _1 is too low in the 
sense that it corresponds to V/a > 10. A choice of V/a ~ 5 
would require an efficiency that is fa 10 times higher than 
their estimate. In addition, there is a factor of fa 3.2 miss- 
ing in the transition from their eq. (15) to eq. (16), since 
(30 kms~ 1 ) 3 /(2G) «3.2Mayr~| / l Mayr" 1 . 

iKrumholz fc Burkertl (|201Ch and IForbes et all (|201lh 
solve the evolution equations for a thin axisymmetric disc 
under the assumption of self-regulated marginal instability. 
The turbulent energy input rate is driven in their model by 
the inflow rate within the disc (corresponding to a fixed £ m 
in our model), and it is dissipated on a crossing timescale. 
They solve the evolution equations for the required inflow 
rate to keep Q — 1, which is very similar to the possibil- 
ity of letting 7^ ch ange in our model (and as assumed in 



6 Note that the factor /J - ^ -2 in cq. (10) of lGenzel et al.l J2008h 
should be corrected to y / /3~ 1 'f~ 2 (Genzel, R., private communi- 
cation). 



ICacciato et al.l 2011), keeping all other fa ctors fixed (Sec- 
tion |43)) . The iKrumholz fc Burkertl (|201Ch solution is very 
similar to ours in case II where we force Q = 1 (Section 
14.111 . However, in their fiducial solution for z fa 2 discs they 
obtain (see their eq. (45)) an inflow rate that has to be 
roughly six times higher than the star-formation rate, or in 
our formalism 7i n f/esF ~ 6. Similarly, the res ulting disc in- 
flow rate at z fa 2 in the IForbes et al.l (|201ll ) fiducial disc, 
where they also assume outflows with 77 = 1, is higher than 
the SFR in disc (Forbes, J., private communication). This is 
another manifestation of our conclusion that the high tur- 
bulence at high redshift cannot be solely driven by the di- 
rect effect of the incoming streams, and it requires addi- 
tional contributions from an intense inflow within the disc, 
and possibly a deplet i on by out flows. We conclude t hat the 
IKrumholz fc Burkertl (|2010h and lForbes et al.l (|201lf ) results 
also agree with the other studies reviewed in this section. 



6 SUMMARY AND CONCLUSIONS 

We have developed an analytic model that describes cer- 
tain aspects of the evolution of galactic gas discs in a cos- 
mological context. Specifically, we addressed the possibil- 
ity that the streaming of external gas into the discs, which 
provides new fuel for disc instability and star formation, 
also directly drives turbulent motions by converting the in- 
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streaming kinetic energy into disc turbulence. The results 
of the model agree with the observations of gas-rich z ~ 2 
galaxies, where a /V TO t ~ 0.2 — 0.3 and Q « 1. In the absence 
of self-regulation, with all the parameters fixed, we find that 
the ratio of turbulent to rotation velocity a/V TO t remains 
constant in time independently of the accretion rate, and 
that discs evolve from violent instability at high redshift to- 
ward stability at low redshift. The constancy of a /V TO t with 
redshift does not agree with the trend suggested by obser- 
vations. In a second case, where the discs are assumed to 
be self-regulated to marginal instability about Q w 1 as a 
result of a duty cycle in the instability and the associated 
star formation and outflows, a /V TO t is found to decline with 
cosmic time in better accord with the observations. In this 
case the duty cycle declines at late times toward long peri- 
ods of stability separating short episodes of instability. Since 
there is no evidence at low redshift for a strong increase in 
the mean timescale for gas consumption into stars, the small 
duty cycle seems to be observationally unfavourable. 

At z ~ 2, we conclude that only if the conversion effi- 
ciency of in-streaming kinetic energy to turbulence is high, 
i.e. close to unity, the direct driving of turbulence by the 
external accretion can be sufficient. Still, unless the SFR ef- 
ficiency is on the high side of the common estimates, or the 
turbulence decay rate is much slower than the dynamical 
time, this scenario relies on disc depletion caused by either 
galactic winds or inflows within the disc, or both. At low 
redshift, the driving of turbulence by the inflowing gas may 
actually be too much if the conversion efficiency remains 
high and the instability and disc depletion are assumed to 
be con tinuously 'on' (see similar conclusions in lDekel et al.l 
|2009bJ). A conversion efficiency or a duty cycle that decline 
with time may help recover the evolutionary trends. 

We wish to emphasize the role played by the depletion 
from the disc in allowing a high velocity dispersion in the 
remaining disc gas. This depletion is a natural result of star 
formation, galactic winds, and inflows within the disc. With 
less gas in the disc, the same energy input would drive a 
higher velocity dispersion. In particular, if stellar feedback 
and/or AGN feedback drive massive outflows, the disc tur- 
bulence will be enhanced even if there is no direct energy 
injection from the feedback source into the remaining ISM. 
If some fraction of the outflows is recycled back into the 
disc, this may add to the direct stirring up of turbulence by 
accretion. 

In a companion paper (|Cacciato et al.ll2011T ). we neglect 
the direct contribution of the in-streaming energy, assuming 
that the disc turbulence is powered by the inflow within 
the disc, which is intimately coupled to the self-regulation 
of the disc instability at Q — 1. We find there that the 
instability is unavoidable at high redshift, because of the 
intense accretion that maintains a high gas fraction, and 
that the discs are driven to Q > 1 at low redshift, primarily 
due to the growing dominance of the stellar component. In 
this model cr/Kot tends to decline at late times. This is in 
qualitative agreement with the trends we find here when we 
impose Q ~ 1. 

In ICacciato et al.l (J201l|) we considered both gas and 
stars in a two-component disc instability analysis, while here 
we limited the analysis to a one-component gas disc, in order 
to allow for an analytic solution. Our results there showed 
that in order to maintain marginal instability, the presence 



of a 'hot' stellar component implies a lower gas velocity dis- 
persion than in the one-component case. This indicates that 
the inclusion of stars in the model considered here would 
have made it even harder to properly suppress the turbu- 
lence driven by external accretion and reproduce marginal 
instability at low redshifts. 

To summarize, our model suggests that turbulence 
driven by cosmological in-streaming may account for the 
high turbulence observed in z ~ 2 discs, but only if the cou- 
pling between this inflow and the disc is high. On the other 
hand, at low redshift our model is in tension with obser- 
vations unless the conversion efficiency of the in-streaming 
energy evolves in a certain way to a low value. Thus, cos- 
mological in-streaming could in principle have an important 
role in driving turbulence in galactic disks, but for this to 
be the primary driver, and to hold throughout cosmic time, 
the energy conversion efficiency between inflow and disk has 
to be fine-tuned. 
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APPENDIX A: CONDITIONS FOR A 
QUASI-STEADY STATE SOLUTION 

For the solution of eq. ([7j| with the assumption Af s i n k = 
M g T _1 to be eqs. (JHJ and @, it is required that M cosmo 
can be treated as a constant, i.e. that it changes slowly with 
respect to the typical timescale of the solution, namely 



dM c , 



i/dt 



™cosmo 

We can write 

dM cosmo /dt _ 

^VJ-COSIXXO 



< l. 



^"cosmo I, 



dJ\L 



cosmo OiVIq, 



dM 



c)t 



dM dt 



(Al) 



(A2) 



since M COS mo changes with cosmic time for a given halo mass, 
while it also changes for a given halo as it builds up its mass. 
For ACDM cosmology we can write approximately 
Afcosmo oc (1 + z) 2 ' 2 oc t~ 1,7 , whe re t is cosmic time 
(JNeistein et al.ll2006l ; [Genel et al.ll2008l ). Therefore, the first 
term equals M c "^a mo 9Af C osmo/<9£ ~ — 1.7t -1 . Further, we can 
replace dM/dt in the secon d term with M COBmo , and define 
t acc = dM cosmo /dM as in iBouche et al.l (|2010h . to obtain 
from eq. (|A1[) the condition 



> \-1.7t' 



t: 



(A3) 



Since the two terms on the right-hand side have op- 
posite signs, considering each of them separately gives the 
most stringent conditions that are required for M cosmo to 
be considered a constant when solving eq. (01. When we 
consider only the growth of M COS mo that accompanies the 
halo growth , eq. (|A3[I becomes t a cc > T. This condition is 
discussed in lBouche et al.l (J2010j), where it is found to hold 
for z < 7. When, instead, we consider only the direct de- 
pendence of Mcosmo on cosmic time, a stronger condition is 
obtained, as eq. (|A3f) becomes i > 1.7r. This is similar to 
the condition i 2> t that is required to obtain the solution 
in eqs. (flO f -ip )! . 

Therefore, we only have left to consider the relation 
between t and r. Since r is shorter than the star-formation 
timescale £sf (as it includes also inflows and outflows), we 
can compare t to £sf as an even more stringent requirement. 
In our model, isF = tdyn/esF, which can be connected to the 
cosmic time if the disc radius is roughly a constant fraction 
A of the virial radius, via the virial dynamical time: fdyn ~ 
O.lAi. For A ~ 0.05, and our fiducial esf = 0.02, we obtain 
t < tsF ~ 0.254. Even if the star-formation timescale does 
not become shorter at higher redshift a s we assume, and is 
instead roughly c onstant isF ~ 2 Gyr (jGenzel et al.l 120101 ; 
iDaddi et al.l 120101 ). r is still short enough compared to the 
cosmic time during the epoch of interest here z ^ 2, or even 
at somewhat higher redshifts. 

In the discussion above, we have referred to the evolu- 
tion of the average accretion rate over an ensemble of halos. 
In principle, the accretion rate of any individual galaxy may 
differ significantly from this average, thus preventing it from 



reaching the steady-state solution. However, both theoreti- 
cal arguments and observational evidence indicate that this 
is likely not the case for the majority of galaxies. From a 
theoretical point of view, an upper limit on the variations 
of the accretion rates in individual galaxies is inferred from 
the scatter around the mean accretion rate found in cosmo- 
logical s imulations, which is w 0.2 — 0. 3dex for either dar k 
matter JGenel et al.l 120081 ) or baryons (|Dekel et al.l 1200981 ). 
Indeed, galaxies in cosmological hydrodynamical simulations 
are often found to follow the steady-state solution for long 
cosmological times (e.g. ICeverino et al.l 120101 ; iGenel et al.l 
120121 : iDekel fc et al.ll2012l ). On the observational side, galax- 
ies populate a rather narrow region in the SFR -stellar mass 
plane with an overall scatter of ~ 0.3dex (e.g. IDaddi et al 



2007J : iNoeske et aT1l2007l : IWhitaker et al.ll2012l ; Tsalmi et al 



20121 ). Such a tight relation could not be explained in the 



presence of large SFR variations. 

We now consider the importance of mergers in driv- 
ing galaxies out of the steady-state solution. It is clear that 
our model does not apply to galaxies after experiencing a 
major merger event, following which they go out of equilib- 
rium in various ways, and in particular may have their disks 
transformed into spheroids. However, major mergers are not 
frequent, as minor mergers and smooth accretio n contribute 
most of the mass in structure formation (e.g. IGenel et al.l 
I2010T ). Minor mergers may trigger the formation of a bar, 
which is however a 'slow' secular process that does not 
affect the processes we consider in this paper to a large 
degree, at least during the violent disk instability phase. 
The y also do not prod uce significant star-formation bursts 
(e.g. Ijoeee et al.ll2009l ). Furthermore, cosmological simula- 
tions show that about half of the external galaxies that 
come in as minor mergers join the disk and share its kine- 
matic s without affecting its steady state (JMandelker fc et al.l 
|2012j). However, even if we consider all mergers with mass 
ratios > 1 : 10 as 'disruptive' as far as our model assump- 
tions go, we should note that galaxies with a stellar mass of 
10 12 Mq (10 11 Mq) at z = have prob ably experienced onl y 
rj 7 (« 1) such mergers since z = 2 (|Hopkins et al.ll2010l ). 
Considering this frequency against the number of galaxy 
dynamical times elapsed within this cosmic time (> 100), 
we can expect most galaxies not to be found in an out-of- 
equilibrium state caused by mergers. 



APPENDIX B: CONVERSION OF POTENTIAL 
ENERGY DURING DISC INFLOWS 

The conversion of potential gravitational energy into turbu- 
lent energy (£ m > 0) is justified by assuming the gas joins 
the outskirts of the disc and generate s turbulence as it mi- 
grates to the bulge throug h the disc JKrumholz fc Burkertl 
I2OIO1 ; ICacciato et al.ll201ll ). rather than by hitting the disc 
from above. In this picture, the inflow rate 7i n f adjusts itself 
to give Q — 1 via a self-regulation loop, since the inflow is 
a result of the instability in the disc. 

In principle, that released energy can be transferred in 
part to the dark matter component, but we do not address 
this explicitly here, as the dark matter is assumed to be 
subdominant inside the disc. The 'immediate' radiation of 
this energy is not considered, because such losses are already 
encapsulated in the dissipation timescale of the turbulence. 
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However, when the energy is transferred to the gas mass 
that stays in the disc, it can also in principle be in the form 
of rotational energy. The purpose of the following analysis is 
to show that under the standard picture of angular momen- 
tum loss, where a small mass gains most of the lost angular 
momentum, the conversion to rotational energy is negligible. 
The angular momentum loss rate is M ln {RV. Let us as- 
sume that angular momentum is transferred into a fraction 
a of the remainder of the disc mass. That fraction of mass 
that is pushed out by obtaining additional angular momen- 
tum has radius R a and velocity V a that can be different 
from that of the disc as a whole. This can be represented by 
writing 

MiniRV = aM g d ( R ^ = aM s RaVa + aMgRaVa (Bl) 
= aM e R a V a + aM g R a V a ^- = \aM % R a V a , 

Ilia I 

where the penultimate equality is based on the equality Va = 
0.5\/GM E i? o /(fi a ) L5 = 0.5V a R a /Ra, which stems from the 
assumption V a = \/GM s /Ra. When the mass aM g gains 
the angular momentum, it is accompanied by an energy gain 
rate that can be written as follows: 



d(GaM 2 /R a 
dt 



GaMl ■ 

~nr Ra 

it IX It a 

2 V 



(B2) 



M int V 2 A 



3V a "Ra' 

This quantity should be compared with the energy that is 
available from the inflows, i.e. M^V 2 . The crucial factor 
is a, which comes in to the final expression in eq. (|B2[) 
via (R/R a ) 2 . If the angular momentum is given to a small 
amount of mass (n< 1) that is expelled to large distances 
due to the angular momentum it acquires, namely R <C R a , 
then the contribution of the inflows to rotational energy in 
the remainder of the disc is small, and most of the energy 
goes to turbulent energy (£ m (a 1). However, if the angular 
momentum of the inflowing gas is distributed evenly across 
the whole disc, then the rotational energy gain that is associ- 
ated with this angular momentum gain is comparable to the 
energy that is released by the inflow, namely little energy is 
left for driving turbulence in the disc and then ( m <l. 



• When gas is accreted at the virial radius, it has poten- 
tial energy of ~ 3V 2 ir relative to the halo centre (for NFW 
haloes), and kinetic energy per unit mass of ~ 0.5V^- We 
assume that most of that poten tial energy is lost t o radia- 
tion (supported by simulations; iDekel et al.ll2009a[ ). Thus, 
at the arrival to the disc, the potential energy between 
the disc edge and the centre is ~ V 2 \ T an d the kinetic en- 
ergy is still w 0.5V;;;. IDekel et all (|2009ah estimate Vj n ~ 
(1.5 — 2)Vvir- We assume in addition V ro t ~ (1 — 1.5)Vvir- 
Thus, u = Vrot/Vin ~ 0.75 ± 0.25. We note that we thus 
assume roughly equal amounts of energy of the incom- 
ing gas, at the arrival to the disc, in kinetic and poten- 
tial form: potential energy of ~ V 2 iv and kinetic energy of 
« 0.5V 2 = 0.5(1.5 - 2) 2 K 2 ir « (1 - 2)Kl. 



APPENDIX C: APPROXIMATIONS FOR 
PARAMETER NUMERICAL VALUES 

The following assumptions are used throughout in order to 
evaluate numerical values for different quantities. 

• The star-formation efficiency per dynamical time is 
small: esf <Si 1. 

• The wind mass-loading factor is not very large: eg F ^> 

(1 + 77). 

• Turbulent energy dissipation is much faster than star 
formation: idis ~ id yn . 

• Angular momentum loss is not much faster than star 
formation: -yi n f ~ esp. Also , one can derive 7j n f S3 0.01 by 
comparing eqs. (19)-(21) in lDekel etHI (|2009bh to the defi- 
nition of Mi n f in Section 12.21 



